Create map of CA counties
Mutate Census Population data for 2020.
%%Mutate Census Population data for 2010. %%```{r} dec_vars_2010 <- listCensusMetadata( name = “2010/dec/pl”, type = “variables” )
smc_pop_2010 <- getCensus( name = “dec/pl”, vintage = 2010, region = "block:*“, regionin =”state:06+county:081“, vars =”group(P1)" ) %>% mutate( block = paste0(state,county,tract,block), pop = P1_001N )
Isolate the blocks for North Fair Oaks 2010.
%%``` {r}
library(devtools)
library(tigris)
smc_blocks_2010 <- blocks("CA", "San Mateo", year = 2010, progress_bar = F)
nfo_boundary <- places("CA", progress_bar = F) %>%
filter(NAME == "North Fair Oaks")
nfo_pop_2010 <- smc_pop_2010 %>%
left_join(smc_blocks_2010 %>% select(block = GEOID10)) %>%
st_as_sf() %>%
st_centroid() %>%
.[nfo_boundary, ] %>%
st_set_geometry(NULL) %>%
left_join(smc_blocks_2010 %>% select(block = GEOID10)) %>%
st_as_sf()
Isolate the blocks for North Fair Oaks 2020.
Place these isolated blocks into a map.
Overlap between 2010 and 2020
%%{r} leaflet() %>% addProviderTiles(providers$CartoDB.Positron) %>% addPolygons( data = nfo_pop_2020, stroke = F, fillOpacity = 0.5 ) %>% addPolygons( data = nfo_pop_2010, color = "red", weight = 0.75, fill = F )